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Abstract. Self-organized critical models are used to describe the 1/f-spectra of rather different physical 
situations like snow avalanches, noise of electric currents, luminosities of stars or topologies of landscapes. 
The prototype of the SOC-models is the sandpile model of Bak, Tang and Wiesenfeld (Phys. Rev. Lett. 
59, (1987) 381). We implement this model on non-periodic graphs where it can become either isotropic or 
anisotropic and compare its properties with the periodic counterpart on the square lattice. 

PACS. 61. 44. Br Quasicrystals - 05.65.+b Self-organized systems 



1 Introduction 

Self-organized criticality (SOC) is observed in various dy- 
namical systems that develop algebraic correlations. The 
systems drive themselves into critical states independent 
of initial states or tuning of parameters. The complex be- 
haviour is characterized by the absence of characteristic 
time or length scales. It is possible to define a set of criti- 
cal scaling exponents that characterize the large scale be- 
haviour and describe the universal 1 /f-power spectra (for 
an overview see Ref. 2 in Q). 

Sandpile automata Q are among the simplest models 
that develope SOC behaviour in the avalanche propaga- 
tion. The model is defined as follows: on each site of a 
lattice, a variable Zi is defined (number of grains of sand). 
At each time step, the variable at a randomly chosen site 
is increased Zi — > Zi + 1. Above a critical value z c , the site 
topples Zi — > Zi — z c and the nearest neighbors j of site i 
are increased Zj — > Zj + 1. Thereby, z c grains of sand have 
been transported from site i to the neighbors j. A top- 
pling can induce nearest neighbor topplings and can cre- 
ate avalanches of arbitrary sizes s. For finite systems, one 
typically uses open boundary conditions (periodic bound- 
ary conditions can lead to never ending avalanches), i.e. 
the sand is just dropping from the table if it reaches the 
border. This introduces the concept of finite size scaling 
(FFS) to compensate cutoff effects. Under the assump- 
tion that FFS is valid the avalanches are described by 
three variables: the number of toppling s, the area a cov- 
ered by an avalanche and the avalanche duration T. The 
probability distributions of these variables are algebraic: 



P(x) = x- T *G{ — ) 

X c 



(1) 



{t x ,Px}- Very recently, the validity of FFS for sandpile 
models was questioned 0j but this discussion is far beyond 
the scope of this article. We refer to || and references 
therein and take FFS as a working hypothesis hereafter. 

The sandpile automata can be isotropic jO, i.e the 
sand is distributed equally between nearest neighbors, or 
anisotropic j|. The question whether both belong to the 
same universality class is still under discussion. Whereas 
real-space normalization group || , field theoretical || and 
late numerical approaches j^j suggest that both belong to 
the same universality class, Ben-Hur and Biham Q ques- 
tion these results. 

In what follows we implement the sandpile model on 
the 8-fold quasiperiodic rhombus-tiling, also called Am- 
mann-Beenker tiling. Due to the specific geometry of qua- 
siperiodic tilings it can become either anisotropic with 
complete non-deterministic toppling rules, or isotropic with 
deterministic topplings, or a mixture of both. 



2 Non-periodic tilings and SOC 



Non-periodic tilings 
try of quasicrystals 



with x £ {s,a,T}. G is a cutoff function and the cutoff 
diverges x c ~ with the system size L — > oo. The 
universality class is determined by the set of exponents 



are used to describe the geome- 
They appear in two distinguished 
classes. The ideal quasiperiodic tilings are produced by a 
cut-and-project-scheme from a higher dimensional (mini- 
mal) embedding latti ce Ec| ,|TT[. Their vertices are given by 
so-called model sets |12| and they show pure Bragg-peak 
diffraction spectra with non-crystallographic symmetries 
(5, 7, 8, ...-fold). Typically, the tilings are built by two or 
more prototiles. Fig. [j] shows the 8-fold Ammann-Been- 
ker tiling consisting of two prototiles, a square and a 7r/4 
rhombus. 

The second class of non-periodic tilings are the so- 
called random tilings. They are constructed by the same 
prototiles as the ideal ones but this time these prototiles 
fit together stochastically in such a way that they fill the 
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entire space, face to face without gaps and overlaps J13| , 
|l4]| . They also can be embedded into a higher dimensional 
lattice but their diffraction properties change. They can 
keep the same symmetry of the ideal counterparts but this 
time, their diffraction spectrum consists of Bragg-peaks 
and diffuse scattering (in 3D) or algebraic peaks (in 2D). 

Beneath the long-range correlations which govern the 
diffraction properties, the main difference of non-periodic 
tilings compared to crystallographic tilings is in their lo- 
cal neighborhood. Whereas we have a fixed translational 
invariant neighborhood for crystals, we have a finite at- 
las of different vertex configurations without translational 
invariance for quasicrystals. For the example of the Am- 
mann-Bccnker tiling, we find 6 different vertices in the 
ideal and 16 different configurations in the random ver- 
sion. The number of nearest neighbors ranges from 3 to 8 
& 

Concerning SOC-models, one now could expect two 
major influences. Both, long-range correlations and the 
diversity of the local neighborhoods of the tilings, could 
develop a difference in the avalanche distribution. It is not 
necessary to test a lot of different quasiperiodic tilings 
because the generic features of quasiperiodicity are the 
same. Symmetry or the type of the tiling plays a minor 
role. Additionally, the local configurations allow different 
implementations of e.g. the sandpile model. For simplicity 
we will focus on three potentially different versions (there 
are more possible but typically they fall in the same cat- 
egories) on the 8-fold Ammann-Beenker tiling: 

(I) anisotropic and non-deterministic: The critical height 
z c and the number of toppling grains n, = z c > are 
equal for all vertices i of the tiling. At every toppling 
step one has to make a random choice to distribute the 
rii grains among the N neighbors (3 < Ni < 8) of the 
toppling vertex i. 
(II) partly isotropic and deterministic, partly not: z c is 
equal for all vertices i but greater or equal to the min- 
imal number of neighbors N m i n — 3 and less than the 
maximal number of neighbors -/V max = 8. The num- 
ber of topplings at site i is n, = min(z c , Ni) where 
Ni again is the number of neighbors of vertex i. This 
way, the toppling rules are locally isotropic and deter- 
ministic for vertices with N < z c but anisotropic and 
non-deterministic for all the others because of the ran- 
dom choice one has to make in order to distribute the 
ni — z c grains among the Ni > z c neighbors. 
(HI) isotropic and deterministic: z c is greater or equal than 
the maximal number of neighbors -/V max = 8, but the 
number of toppling grains of vertex i is always n, = Ni. 
This way, every bond will transport the same amount 
of grain in both directions and the number of topplings 
is given locally by the number of neighbors. 

After a suitable thermalization, we perform 2 * 10 6 
(L = 16, 32, 64) and 5 * 10 5 (L = 128, 256) avalanches 
for the different scenarios. Each avalanche is seperated 
by an appropriate thermalization time to prevent cross- 
correlations. The system size (L x L) is given by the length 
L which varies from 16 to 256 whereas the absolute scale is 




Fig. 1. The ideal Ammann-Beenker tiling and a critical state 
at the end of an avalanche of type (I). The symbols on top of 
the vertices code the number of grains of sand z, and z c = 5. 
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Fig. 2. Comparsion of anisotropic avalanche distributions of 
an ideal and a random tiling version of the Ammann-Beenker 
tiling at L=128. 



given by the bond length of the tiling that is equal to one. 
For each setting, the avalanche distributions are stored. 
According to the finite size scaling assumption, the data 
for the different length scales should collaps to a single 
curve that is determined by the critical exponents t x , [3 x . 
Fig. H shows the data collaps of the probability distribu- 
tions of the number of topplings s for version (I) , Fig. ^ 
for version (II), and Fig. | for the isotropic deterministic 
version (III). Table | gives the extracted numerical values 
for the exponents. The errors given in the brackets are the 
statictical errors of the fitting procedure. 

Within the limits set by the small statistic, all three 
version belong to the same universality class: the crtiti- 
cal exponents and also the scaling functions seem to be 
the same. If one compares the distributions of the ideal 
quasiperiodic tilings and the corresponding random tiling, 
there is no apparent difference (see Fig. ||). Looking at the 
crystalline counterpart, the 2D square lattice, the state 
of the art values for the exponents are t s = 1.27(1) and 
(3 S — 2.73(2) These values are determined by massive 
parallelized computation. From the difference to Table |l| 
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Fig. 3. Data Collaps of the probability distribution of the 
number of topplings s for the anisotropic, non-deterministic 
version (I) on the ideal Ammann-Beenker tiling. System 
lengths are L=16 (o), 32 (□), 64 (o), 128 (A), and 256 (+). 
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Fig. 4. Data Collaps of the probability distribution of the 
number of topplings s for version (II) on the ideal Ammann- 
Beenker tiling. System lengths and symbols as in Fig. 3. 



10"- 




10 

s/L J 



Fig. 5. Data Collaps of the probability distribution of the 
number of topplings s for the isotropic, deterministic version 
(III) on the ideal Ammann-Beenker tiling. System lengths and 
symbols as in Fig. 3. 



Table 1. The citical exponents t 3 ,/3 s of the three tested ver- 
sions (Statistical errors are given in brackets). 



Version (I) 


Version (II) 


Version (III) 


t s 1.20(3) 


1.19(3) 


1.19(3) 


Ps 2.5(2) 


2.5(2) 


2.5(2) 



one cannot conclude that there exists a different univer- 
sality class for quasiperiodic graphs. The values of Table 
[j] are more comparable to older estimates for the square 
lattice based on similar statistics The distributions 
for the area a and duration T are not shown because they 
behave in the same way. 



3 Conclusion 

We have performed numerical estimates of the critical ex- 
ponents of the 2D sandpile model on the 8-fold quasiperi- 
odic Ammann-Beenker tiling. Thereby, three different ver- 
sion (anisotropic-non-deterministic, isotropic-deterministic 
and a mixed one) were implemented on the ideal and on 
the random Ammann-Beenker tiling. The measured expo- 
nents are comparable to early numerical results of the 2D 
square lattice and might suggest that all tested version 
belong to the same universality class as the 2D square lat- 
tice. The problem however remains to substantiate this, 
e.g. by better statistics. This can only be achieved by a 
massive use of a parallel computing which is not at hand 
at the moment. But the situation seems to be similair 
to the investigation of Ising-models JlTj which led to the 
conlcusion that aperiodic and periodic versions belong to 
the same universality class. 

The author wants to thank S. Maslov and M. Baake 
for discussion and helpful comments. 
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